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a^ : 

, I In this paper, we study generic non-local NJL models for four-dimensional QCD 

' in the limit of a large number of colours. We diagonalise, through a Bogoliubov-type 

T— I ■ 

^__{ ■ transformation, the Hamiltonian of the model completely in terms of the compound 

o^ : 

, operators creating/annihilating mesons — bound states of dressed quarks and an- 

I tiquarks — and demonstrate a one-to-one correspondence between the Bogoliubov 

p • diagonalisation condition and the Bethe-Salpeter equation for bound states. We 

Oh. conclude that, in the leading order in Nc, the nontrivial contents of the theory is 



X 



entirely encoded in the chiral angle — the solution to the mass-gap equation for 
dressed quarks. As a consequence, we extend the statement of existence of excited 



■ solutions to the mass-gap equation, called vacuum replicas, beyond the BCS level 

and prescribe a new index to the mesonic operators — the label of the vacuum state 
in which the corresponding meson is created. 

PACS numbers: 12.38.Aw, 12.39.Ki, 12.39.Pn 



I. INTRODUCTION 

For nearly half a century Nambu-Jona-Lasinio (NJL) |l| class of models 0, 0, 0] have 
been quite useful tools for the on going study and understanding of low-energy phenomena 
in QCD. Quark models with non-local quark current-current microscopic interactions, pa- 
rameterised by kernels of particular forms, constitute prominent steps towards building an 



2 



effective low-energy theory for the strong interactions. Besides regularising the ultraviolet 
divergences of the theory, these nonlocal quark kernels bring in the necessary interaction 
scale needed to make contact with the hadronic phenomenology. Models with quark instan- 
taneous interaction modeled after the harmonic oscillator kernels l2Ly,y| or kernels obtained 
in the Coulomb-gauge QCD as well as a variety of other approaches based on the lin- 
ear confinement [6|, constitute as many examples of closely related approaches exploiting 
the idea of chirally nonsymmetric solutions of a nonlinear gap equation. In this paper, we 
continue our studies of quark models ^ with instantaneous inter-quark interaction. An im- 
portant feature of this class of models is the essential covariance of results and conclusions 
across a variety of possible forms for the confining kernel. For example, the pion-pion low- 
energy scattering problem can be solved in the general form, in terms of the pion mass and 
decay constant, regardless of the explicit form of the quark kernel, provided it is confining 
0]. Another example of such universality is given by the existence of multiple chirally non- 
invariant solutions to the mass-gap equation, predicted in 0, 1^, found both analytically 
and numerically for all allowed power-like confining potentials in 1^, and independently 



confirmed by calculations done in a different approach in Ul!|. The same conclusion of uni- 
versality holds for any low-energy phenomena considered in the framework of such class of 
models and, therefore, we dare to expect to be able to reproduce the low-energy properties 



of QCD. An important progress was achieved in 1^ where a two-dimensional model for 
QCD 3] in the axial gauge was considered and a second, Bogoliubov-type, transforma- 
tion was performed in the mesonic sector to fully diagonalise, up to corrections of order 
0{1/ \/Nc)i the Hamiltonian of the theory in terms of the mesonic creation/annihilation 
operators. The Hamiltonian approach to two-dimensional QCD in the light-cone gauge was 
also developed in In the present paper, we extend this result to four-dimensional QCD 
and establish a one-to-one correspondence between the Bogoliubov condition (which ensures 
that the Hamiltonian is diagonalised in the mesonic sector) and the Bethe-Salpeter equa- 
tion for bound quark-antiquark states. We build the corresponding amplitudes explicitly 
and study in detail the cases of the chiral pion and of the p-meson. We conclude that 
solving the Bethe-Salpeter equation in the ladder (rainbow) approximation is equivalent, in 
the leading order in Nc-, to go beyond BCS level and to diagonalise the quartic term in the 
Hamiltonian of the theory. This feature is general and holds for any confining quark kernel. 
We use the oscillator-type quark kernel just as an example to illustrate the method. We 
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also conclude that all the information needed for the second and final mesonic Bogoliubov 
transformation, is already contained in the solutions of the mass-gap equation appearing at 
the BCS level. Therefore, in the leading order in Nc-, the chiral angle — the solution to the 
mass-gap equation — remains the only nontrivial characteristic of the theory and it defines 
the latter completely. As a result, we discuss the multiple chirally noninvariant solutions of 



mass-gap equations — the vacuum replicas 



Id — and conclude that, once such solu- 



tions exist, this information gets carried beyond BCS and appears as an extra spectroscopic 
index for mesonic states on top of the usual set {P, n, J^^} — namely, each meson carries 
the following quantum numbers: the total momentum P, the radial excitation number n, 
the J^^ set, and, finally, the label of the vacuum state in which it is created. We discuss 
possible consequences of this generalisation. 

The paper is organised as follows. The second section is first devoted to the diagonal- 
isation of the Hamiltonian in the quark sector (subsection A) and then to its subsequent 
diagonalisation in the mesonic sector (subsection B). We end this chapter with the discussion 
of the equivalence between the mesonic Bogoliubov- Valatin transformation and the Bethe- 
Salpeter equation. Vacuum replicas are studied in the third chapter. We give the necessary 
details of the theory of replicas, in the first subsection, and comment on the tachyon prob- 
lem, in the second subsection. Finally, in the third subsection, we generalise the theory of 
replicas beyond BCS and discuss the hadronic spectrum in the replica vacua. The fourth 
and last chapter of the paper is devoted to our conclusions and outlook. 

II. DIAGONALISATION OF THE QUARK HAMILTONIAN 

A. Diagonalisation in the single-quark sector and the mass-gap equation 

In this section, we give the necessary details of the quark model P, 0, 1^ to be used in 
this paper and describe the Hamiltonian diagonalisation in the single-quark sector (BCS 
level). Naively one should expect this to be the end of the story. However, as we shall 
see, the requirement of confinement, that is, no free asymptotic quarks may be allowed, will 
completely change this picture and we shall end up with a diagonal Hamiltonian written in 
terms of asymptotic mesons rather than asymptotic quarks. How this comes about will be 
the subject of this paper. 
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The Hamiltonian of the model, 

H = j d'xijix,t)(^-t^-^ + m}ji;ix,t) + ^ j di'xd'y J^{x,t)K;l{x - y)Jtiy,t), (1) 

contains the interaction of the quark currents J^{x,t) = ^jj(x,t)'j^^ip{x,t) parameterised 
through the instantaneous quark kernel, 

Kf^i^-y) = g^og^oS^'Voilx-yl)- (2) 

In Eq. ((21) we have used the simplest form of the kernel compatible with the requirement of 
confinement. Spatial components can be also included in the kernel (j2I) fl- They affect, for 
example, the form of the mass-gap equation for the chiral angle and spin-dependent terms 
in the effective quark-quark interaction after a Foldy-Wounthuysen transformation (see, 
for example, papers 15| where a heavy-light quarkonium was considered in the modified 
Fock-Schwinger gauge and an effective inter-quark interaction was derived). 

But it is a fact, already mentioned in the introduction, that despite all the aforementioned 
detail variations of possible quark kernels, they all share with the kernel (j21) the same low 
energy properties (Goldstone pion, vr — vr Weinberg scattering lengths, Gell-Mann-Oakes- 
Renner formula, and so on) so that it is sufficient to stick to the form Q as the mainstay 
of this work and comment on the generalisations when it is the case. Finally it is suffi- 
cient for phenomeno logical applications that Vo(|a;|) should interpolate between long range 
confinement and short-range Coulomb behaviour. An example of such a phenomenological 
kernel which is able to address both, heavy-quark and light-quark limits of the theory, was 
considered in detail in ^. A pure confining kernel of the generalised power-like form, 

Vo{\x\) = KS+'\xr, (3) 

with Kq being the mass parameter, was studied for the first time, in the series of papers P] 
and an analysis was performed concerning the general properties of the mass-gap equations 
for such potentials. This analysis was extended in a recent paper [lo| and the range of 
possible powers a was identified to be ^ a ^ 2. In what follows, we do not need to specify 
any particular form for Vo(|x|) other than requiring it to be confining. When needed, we 
shall use the oscillator- type quark kernel corresponding to the case of a = 2 in Eq. Q, 
to exemplify some general results. Convenience of such a choice follows from the fact that 
the Fourier transform of the potential x"^ is given by the Laplacian of the three-dimensional 
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delta function, so that all integral equations can be transformed into differential equations 
which are much easier for analytical and numerical studies Q, 0, Oj ■ 

Now we proceed in the standard way by defining dressed quarks with the help of a 
Bogoliubov-Valatin transformation having the chiral angle if{p) = (pp as a parameter: 




1 

1 
V2 



m(0), 



(5) 



y/l + sin(y9p + ^1 - sinv9p (ap) 
1^1 + siny^p - yT^^sm^ (a^ 

t) = e'^-'Up, 0), t) = e^-'d^i^-p, 0), (6) 

where Ep stands for the dispersive law of the dressed quarks; a being the colour index, 



a = l,Nc- It is convenient to define the chiral angle varying in the range — ^ < ^ | 
with the boundary conditions 93(0) = |, (p{p 00) 0. Notice that this definition is not 
unique and in some cases one is forced to choose ip{0) = — |, as we shall see. 
The normal ordered Hamiltonian can be split into three parts, 

H = -E'vac+ : H2 : + : H4 '■, (7) 

and the usual procedure is to demand the quadratic part : H2 : to be diagonal or, equiva- 
lently, that the vacuum energy E^^c should become a minimum. Then, the corresponding 
mass-gap equation ensures the anomalous Bogoliubov terms wd^ — db to be absent in : ^ 

Ap cos ipp — Bp sin ipp = 0, (8) 

where 



1 f d'^k 

Ap = m + -Cf J j^-^Vo{p- k)smipk, (9) 
1 f d^k 

Bp = p + -Cf J {pk)Voip- k) costpk, (10) 

Cf = ^{Nc—l/Nc) being the SU {Nc)c Casimir operator in the fundamental representation. 
The large- A^^c limit implies that the product CfVq remains finite as Nc 00, so that, for 
the power-like form (jHl), an appropriate rescaling of the potential strength Kq is understood. 
In the remainder of this paper we shall absorb the coefficient Cf into the definition of the 
potential, CirVo(|x|) = by a trivial redefinition of Kq. 
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As soon as the mass-gap equation is solved and a nontrivial chiral angle is found, the 
Hamiltonian ((Tj) takes a diagonal form, 

H = E^^, + j2Y. [ 4^,Ep[hiMKs{p) + du-p)du-p)i (11) 

and the contribution of the : : part is suppressed as 1/y/Nc- The dressed quark dispersive 
law becomes 

Ep = Ap sin ipp + Bp cos ipp, (12) 

and this completes the diagonalisation of the Hamiltonian in the quark sector to order 
l/VN^ (the BCS level). 

For further convenience we write out both the mass-gap equation and^he dressed quark 
dispersive law Ep for the harmonic oscillator potential, V^(|x|) = Kgl^P 







sin ipp = [pVp + 2pv5p + sin 2ipp] + mp^ cos ipp, (13) 



Ep = m sin (Pp + p cos (pp — Kq 



V^p ^ cos^ ipp 



(14) 



2 ' p2 

In the remainder of this subsection, if not stated otherwise, we consider the chiral limit 

171 = 0. 

We end this introductory section by stating two general properties held by mass-gap equa- 
tions, Eq. (jHl). First of all, notice that using the explicit form of the dressed quark dispersive 
law Ep (fT^ . the mass-gap equation can be rewritten in the form of a Schrodinger-like 
equation with the zero eigenvalue, 

[-KlAp + 2Ep]tP = 0, (15) 

where if) = sin (pp. The generic mass-gap equation (jH} for an arbitrary power-like confining 
potential Q admits the form similar to ()15|) [lo| . 

The second general property concerns the asymptotic behaviour of the solutions of the 
mass gap equation (the chiral angle) at large momenta. Taking the general form of the 
mass-gap equation for the power-like potential (jHl) [lfll |. 

3 . T^a+irf , IN • ™ /"°° ( pk sm[<pk - <pp] cos <pk sin <pp \ 

p sm^, = K,- Via + 1) sm - j ^ - | |^ _ ,,,,, + _ _ | ' (16) 
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and performing an expansion for p — ^ oo, we get ^ 

From the behaviour p7|) we see that, in the chiral hmit, the sign of the solution Lpp at 
large momenta is opposite to the sign of the chiral condensate. Beyond the chiral limit the 
expression ()17|1 defines how the chiral angle approaches the free large-momentum asymp- 
tote arctan^. From the Gell-Mann-Oakes-Renner formula |l6| one can relate the chiral 
condensate with the pion mass M^-, 

flMl = -2m{qq), (18) 

so that, according to Eq. (fTTj) . becomes positive if the chiral angle approaches its large- 
momentum asymptote from above, and negative if it approaches from below. In the latter 
case the pion becomes a tachyon and this situation requires a special treatment which will 
be discussed in detail in chapters II and III. 



B. Diagonalisation in the mesonic sector and the bound-state equation 

1. Quark Hamiltonian in the space of quark- antiquark pairs 

In this subsection, we go beyond BCS level and include the quartic part : : of the 
Hamiltonian (0) into consideration. 

After diagonalisation at BCS level, we could have used the formal invariance, under 
arbitrary Bogoliubov-Valatin transformations (BV), of the quark field ipa{Xit) — here un- 
derstood as an inner product between a Hilbert space, spanned by the spinors Us{'p) and 
Vsi—p), and a Fock space, spanned by the quark annihilation and creation operators. This 
invariance requires that any given BV rotation in the Fock space, must engender a corre- 
sponding counter-rotation in the associated spinorial Hilbert space so as to maintain the 

^The case of the harmonic oscihator potential, a = 2, has to be considered separately. As immediately 
follows from Eq. (|15(l and the free limit of the quark dispersive law, Ep « p, the asymptotic behaviour of 

p — *oo 

the chiral angle in this case is given by the Airy function. Still the qualitative conclusion made below holds 
true. 
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inner product ipa{x,t) invariant. Thence, the new BV-rotated spinors will carry the infor- 
mation on the chiral angle and, therefore, so does the quark propagators. Equipped with 
these new Feynman rules, we can then proceed to evaluate a wealth of physical processes 
among which the Bethe-Salpeter equations for bound states play a central role. Formally, 
the Hamiltonian of Eq. (^, when written in terms of the quark fields ipaixjt), is invari- 
ant under such BV transformations and so it will contain — through its : : term — 
Bogoliubov anomalous terms of various types. 

It will turn out that the set of all these Bethe-Salpeter equations, each one of them as- 
sociated with a corresponding mesonic bound state, with their right mixture of positive with 
negative- energy amplitudes, are exactly what is required in order to get rid all the Bogoliubov 
anomalous terms with just four quark creation(annihilation) operators. The remaining Bo- 
goliubov terms (responsible for generic hadronic decays like A ^ B + C ) will be suppressed 
as 1/y/Nc- This scenario is quite plausible once we realise that a group of four-quark cre- 
ation(annihilation) operators — two quarks and two antiquarks — corresponds to a bilinear 
creation (annihilation) of a group of two mesons and, as we know, these bilinear anomalous 
terms can be made to disappear by the use of an appropriated bosonic Bogoliubov- Valatin 
canonical transformation. 

For the implementation of the diagonalisation program we shall follow the method of 
refs. I12I] generalising it to the four-dimensional case. The basic idea of the method is the 
following. At BCS level, we consider quarks as though they were free, with confinement 
to be implemented at later stages, when building the Bethe-Salpeter equation for bound 
states. Instead we now start with the ab-initio confinement requirement that only quark- 
antiquark pairs are allowed to be created and annihilated, not single quarks, and that the full 
Hamiltonian ((7j) must take a diagonal form in terms of the operators creating/ annihilating 
whole quark-antiquark mesons. Therefore, every time a quark (antiquark) is created or 
annihilated, an accompanying antiquark (quark) must be created/annihilated as well. To 
write this ab-initio condition, let us introduce a set of four colourless operators: two operators 
which count the number of quarks and the number of antiquarks. 



BsAp,f) = ^j:JUp)kAif), D,Am = ^j2ji(-p)i,(-ji), (19) 



MUp,if) = ^j:JUi^)dU-p), MsAm = ^j:j^si-p)kAf)- (20) 





and two operators which create and annihilate a quark-antiquark pair. 
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The coefficients in the definitions ()19|) and ()2U|) are fixed in such a way that, in the large- A'^c 
limit, all four operators obey the standard bosonic algebra, with the only nonzero commu- 
tator being 

[MsAv.f) MlAl Q)] = {27rf6^'\p- m^fS^'Hp - q)SsJs'a'. (21) 
The pair (|T9|l is sufficient to express the leading part of the Hamiltonian (|TT|l . 

H = E,,, + J2 J -0j5^piBssip,p} + Dssip,p}l (22) 

whereas the suppressed part : H4 : requires all four operators. 

According to our ab-initio confinement requirement, the operators B and D cannot be 
independent and have to be expressed through the operators and M. One easily builds 
such relations to be jlJ] 

73^// 



\/Nc ^ J (2vr) 



'3.// (23) 



The anzatz ()23p is the most crucial point of the method. Formally, one can also view it 
as another solution to the set of equations given by the commutators of the operators (fTTHl 
and (j2(J|) . Notice that the relations (j23j) correspond to the "minimal substitution", that is, 
each quark is accompanied by only one antiquark, and vice versa, whereas, in the general 
case, a whole quark-antiquark cloud should be also created together with the accompanying 
particle. Since the creation of each next quark-antiquark pair in the cloud is suppressed 
by an extra power of l/Nc, then Eq. represents just the leading terms of the full 
expressions. 

After substitution of the anzatz (j23|) into the Hamiltonian (|7j) we have no longer a general 
suppression of the quartic part of the Hamiltonian as compared to the quadratic part. Only 
a number of terms in : : remain suppressed and will be omitted henceforth. The resulting 
Hamiltonian, in the space of colourless quark-antiquark pairs, takes the form: 

^ = ^-^ + / 1^^^^)' 
where we have separated out the quark-antiquark cloud centre-of-mass motion. 
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For the sake of simplicity, let us consider the Hamiltonian density Ti at rest (P = 0), 
n^nP = 0) = J2[ ^,2E,MUmMs2sMp) + ^ E / ^s^s^P- ^ 

+ [v'^{p,q)\s^szs^S2Msrs2{P,P)MsiSi{q,q) + {P, Q)\s^S3SaS2^S3sAP^P)^I^S2'^^^ ' 

with amplitudes v given by 

K+(p,g)]sis3S4S2 = [u\M'^sM]HS-^^s2{.-p)l 

[^"+(p,g)]sis3S4S2 = [vI^{-v)usM]HS-^Us2{:p)1 
[v—{p,q)]s^s^s^s2 = H,i-Q)vs-A-p)]HM^^2iQ)]- 

Although there are only two independent amplitudes in ()26|) . for example, v~^~^ and , 
with the two remaining amplitudes, v and v~^, easily related to the first two by Hermitian 
conjugation and renaming of spin indices, we prefer to keep all four amplitudes so as to give 
a more symmetric form to both bound-state equations and the effective diagrammatic rules 
to be derived shortly. 

Notice that in general, the vacuum energy in the Hamiltonian ()24|) . E'^^^, differs from the 
vacuum energy i^^ac in the Hamiltonian (^TJ) by terms coming from the commutators of the 
operators M and M'^ . 



2. The case of the chiral pion 

Before we proceed with the diagonalisation of the Hamiltonian ()25j) in general form, let us 
consider the case of the chiral pion as a paradigm for the general case without unnecessary 
complications to cloud physics. The spin and spatial structure of the pion — a. ^Sq state 
with J = L = S = — can be chosen in matrix form so that the operator Mss/{p,p), 
creating the pion at rest, is represented as 

Substitution of Eq. (j?fj) into the Hamiltonian gives a simpler expression in terms of the 
radial operators M{p) and M^{p), 



(27) 
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+T+-{p, q)M\q)M\p) + T~+{p, q)M{p)M{q) + T^' (p, q)M\q)M{p)j , (28) 

where the amphtudes T^^ {p, q) are trivially built from the amplitudes v^^ {p, q) (see 
Eq. ()2()|l ). after inclusion of the potential V{p — q), integration over angular variables, and 
taking a spin trace against the pion spin wave function k, as it was defined in Eq. ()27j) . 

To have a fully diagonalised Hamiltonian it is now sufficient to perform a Bogoliubov 
transformation, 



I Hp) -- 

\ M^p) 



(29) 



or, inversely. 



ml 



p^dp 

(2^ 
p^dp 



M{p)vt{p)-M\p)v-{p) 
M\p)vt{p)-M{p)^-{p) 



(30) 



(27r)3 

where the operators m\ and m^r create and annihilate pions at rest. 
Using the commutator of the radial operators, 

[M(p), M\q)]=^^8{p-q), (31) 

following from Eq. (PT|) . we find that 

[r^., rht] = I ^ [^+2(p) - ^;2(p)] . (32) 

Therefore, requiring the canonical commutator between the pion creation and annihilation 
operators to be [m7r,mj.] = 1, ensures that the Bogoliubov amplitudes ^9^, playing the role 
of the pion wave functions, must in turn be subject to the normalisation ^ 



^Although Eqs. seem not to follow from Eqs. H29(l directly, which is a consequence of the naive 

truncation we have made leaving only the pion instead of the whole tower of mesonic states, in fact, Eqs. H29|l 

and H3U|I are self-consistent with each other at the given level of truncation. See Eqs. H64|l and Ht)5|l below. 
■^Notice that the wave functions can be chosen real. This holds for all observable mesons as well. If the 

mesonic wave function still contains an imaginary part then Eqs. H29|) - (|33|l should be changed accordingly, 

so that, for example, ip^"^ in the normalisation condition H33|l should be changed for |<y5^p. Below we discuss 

the situation with the tachyon where one is forced to deal with an imaginary mass and, as a consequence, 

with its imaginary wave functions. 
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(34) 



To complete the program of bosonic diagonalisation of the Hamiltonian (0), we need 
only to look for a pair of ^p^{p) to diagonalise the hadronic Fock subspace spanned by zero 
momentum pions, |vrp=o) = "^tl^); 1^) being the mesonic vacuum (see the discussion below). 

Now, if ftip) v^~(p) are made to obey the following set of equations, 

[2E,-M^]^t{p) = J ^,[Tt^{P.q)^t{<l) + Tr{p,q)^-{q)] 

[2E, + M^]^-{p) = I ^,[T;-'{p,qWM +T"{p,qW-m 

we can, by a simple substitution in Eq. (j2Hl), get rid of the anomalous Bogoliubov terms to 
have (n|7i^|7rp=o7rp=o) = and {ttp=oHp=q\'Ht,\Q) = 0. The non-anomalous terms can be 
readily summed to yield: 

= M^mlm^ + . . . , (7rp=o|^7r|vrp=o) = M^, (35) 

where the ellipsis denotes terms suppressed by Nc- For example, the quartic term in the 
pion creation/annihilation operators responsible for the pion-pion scattering can be identified 
among the latter. 

Notice the difference between the quark BCS vacuum |0), annihilated by the dressed quark 
operators b and d, and the mesonic vacuum annihilated by the mesonic operators, for 
example, by 171^^. The two vacua are related by a unitary transformation, |0) = with 
the operator W creating pairs of mesons and taking such a form that, for example, 

m^\n) = m^U^\0) = U^{Um^U^)\0) oc U^M{p)\0) = 0. 

Since creation of any extra quark- antiquark pair is suppressed in the large- A^^c limit, so is 
the deviation of the operator f/"'' from unity. As a result, the chiral condensate calculated 
at BCS level, using the BCS vacuum |0), coincides, up to I/Nq corrections, with the full 
quark condensate calculated in the mesonic vacuum 

It is easy to recognise in Eq. the Bethe-Salpeter equation for the pion studied in a 
number papers 0, Q| . To see this equivalence we proceed in five steps. 
Step 1: 

First let us define dressed quarks and derive the mass-gap equation. To this end we 
consider the quark mass operator in the rainbow approximation, 

zEip) = J ^,V{p- k)TSp{k, ko)T, (36) 
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where SFip,Po), 



Sf{p,Po) 



A+(p)7o ^ A-(p)7o 



Po- Ep + ie po + Ep-ie' 



(37) 



(p) = ^ [1 ± 7o sin ipp ± (dp) cos v^p] , 



is the Feynman propagator of the dressed quark with the dispersive law Ep. F is a generic 
Dirac matrix defining the Lorentz nature of the confining interaction. Parameterising T,{p) 
and Ep as 

= [Ap - m] + (7p) [Bp -p], Ep = Ap sin ipp + Bp cos ipp, (38) 

we find that 

Sp^iPo^P) = loPo - {.lP)Bp - Ap. (39) 
Taking the integral in the energy in Eq. (jHBj) . we arrive at the self-consistency condition, 

Ap cos ipp — Bp sin = 0, (40) 

which is the mass-gap equation for the chiral angle (pp. For F = 70 the auxiliary functions 
Ap and Bp are given by Eqs. and (fTnjl . Notice, however, that the general form of Eq. pUjl 
holds for any Lorentz nature of confinement (see ^ for a detailed discussion). This completes 
the matching, at BCS level, between the Hamiltonian and diagrammatic approaches to the 
theory. 
Step 2: 

Proceed beyond BCS level and consider the Salpeter equation (in the Dirac representa- 
tion) for a generic meson in the rest frame, 



X{p; M) 



-Vip- q) VSp^q, go + M/2)x(g; M)^^(g, go - M/2)F, 



(41) 



where xiP'i stands for the mesonic Salpeter amplitude. 



Step 3: 

Due to the instantaneous nature of V{p — q), it is a simple matter of poles logistics to 
see that, upon integration in the energy, the only surviving combinations of the Feynman 
projectors A^*" entering in Eq. (jH} are (A+7o)x(A^7o) and (A^7o)x(A+7o), so that we can 
decompose x i^i two distinct amplitudes, x'^^ aiid To this end we get rid of the energy 
denominators in Eq. (j^T| by taking the integral in the energy. 



dqo 

27Ti 



qo±M/2- Eq + ie 



1 



qo T M/2 + Eq-te 



2E„tM' 



(42) 
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and define: 



V W in- M) = vH (q- M) = 

Then the Bethe-Salpeter equation (PT|) amounts to a system of two coupled equations: 



[2E, - M]xW = - I ^^(P - r [(A+7o)x'+kA"7o) + (A-7o)xf^(A^7o)] T 



(43) 



Step 4: 

Sandwich both equations in between the quark spinors and use the definition of the 
projectors, 

A+(p) = 5^^isi(p)®<(p), A-(p) = ^t;,,(-p)®t;t^(-p), (44) 

S1S2 SiS2 

to cast Eq. (jl3|l as 



S3S4 



[2Ep + M] [Vs,X^~^Us^] = J -^;^Vip- Q) {[Vsr'i^Us3][Us3X^^^Vs^][VsJUs^] 



Step 5: 

Define $+,2 = [usiX^^^Vs^] and, similarly, ^-^^ = [vs.x^'^Us^] to get: 

[2E, - Af]<|.+ ,^ = - E / tC^^^^"- {y''Us,s,s,Ks, + [v'-Us.s^s.KsA 

S3S4 ^ ' 



(45) 



(46) 

In the coefficients f one can easily recognise the amplitudes which appeared in the 
Hamiltonian (j25|) . For the case F = 70 they are given in Eq. ()26p. Further details and 
explicit forms for the amplitudes ()26|) in terms of the chiral angle can be found in refs. . 
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For establishing graphical rules it is convenient to include the potential into the definition 
of the amplitudes and to rewrite them in the form: 



S1S2S3S4 
[T-^{p,q)] S1S2S3S4, 

[T''{P,Q)Us2S3Si 



[msi {p)Tvs3{-q)] [-V{p- q)] [Us4 {qjTvs., {-p}] , 
[Vs, ( -p) Tus, {q)][-V{p- q)] [vs, {-q)Tvs 



(47) 



-Q)'^vs3i-P)][-Vip - q)][Us^{p)Tus2iq)], 



or simply, 



(48) 



T++ = [uTu][-V][vTv], T+- = [uTv][-V][uTv], 
T'+ = [vTu] [-V] [vTu] , T— = [vTv] [-V] [uTu] . 

In Fig. 1 we depict the steps leading to Eq. (j46|) with the amplitudes (jUj). In this figure, 
the first and the second lines represent the r.h.s. of the first equation in the system (jlHI), 
sandwiched between the quark spinors u and f , and the r.h.s. of the first equation in the 
system fl46|) . respectively. Similar diagrams can be drawn for the second equation of both 
systems — the only difference will be the ordering of the spinors u and v. Thus we introduce 
a diagrammatic technique relating Bethe-Salpeter amplitudes and vertices to the coefficients 
of the Hamiltonian in the representation of mesonic operators. 

Now we need to project the spin-angular wave functions ^^-^s^iP) onto the wave function 
with the definite total momentum, parity, and charge conjugation number. For the pion this 
is quite easily done. We have 

1 



V2 



0-2 



Foo(p)^^(p), Y,o{p) 



S\S2 



(49) 



which should be compared with Eq. (j27j) . Finally, we need only to perform a spin trace and 
to introduce the T^^ amplitudes for the pion, such as, 

T++(p,g) = -^ J dn,dn,Y*,iff)Vip- q)YooCq)Sp{a2[m^um'^2[vi-q^^^^ (50) 

to arrive at the pion Bethe-Salpeter equation in the spin representation, which coincides 
with the bosonic mass-gap equation (jS^)- For the case of the harmonic oscillator potential 
it can be rewritten in a matrix form, 



1 






+ ^0 


1 



COS^ ipp 

2 p2 





1 1 




1 










1 1 




-1 






Kip) 



= 0, 
(51) 
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FIG. 1: Going from the Dirac representation to the spin representation: one uses the definitions of 
the projectors A='=, Eq. H44() . (dash-dotted hne) to reshuffle the spinors towards the F vertices (soUd 
line). The Salpeter amplitudes x'^^ are vectorially multiplied by the adjacent spinors in order to 
construct two functions — two-by-two matrices in spins: and (fat solid line). 

for the radial wave functions i^^{p) = P'^tip) which are rescaled so as to obey the one- 
dimensional normalisation ^, 

J dp[uf{p)-u-\p)] = l. (52) 

Comparison of the Hamiltonian (j^Hj) with the bound-state equation pUj) allows us to estab- 
lish a dictionary between the terms in (PHjl and the diagrammatic rules used in the literature 
in order to arrive at the Bethe-Salpeter equation (jl^jl . These rules are easily deduced from 
Fig. 1 and will be exploited in what follows in order to diagonalise the full Hamiltonian ()25j) . 

Notice several important properties of the bound-state equation ()34|) and its solutions. In 
the chiral limit the pion mass vanishes, = 0, and the pion wave functions are real and are 
related to one another and to the chiral angle as (p^{p) = — (/3~(p) = AC simpp, AC being the 
norm. Let us use the simple case of the harmonic oscillator kernel to illustrate this. Using 
the explicit form of the T-amplitudes taken from the papers [J], one easily finds that, for the 
case of the state, T++{p, q) - T+-{p, q) = T-+{p, q) - T--{p, q) = ^^6{p - q)Kl\, 
where the Laplacian is actually reduced to its radial part (see Eq. (jSH)). Substituting this 
result in the bound-state equation (pij) for v^J(p) = — V5~(p) = ^-k, and using ik^ = 0, we 
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get 

[-iro3Ap + 2Ej^^ = 0, (53) 

which coincides with the mass-gap equation p5|) with ip^^ = ip = siny^p. Derivation of the 
pion bound-state equation as a Shrodinger-hke equation for the generic form of the potential 
is given in Appendix A. 

Thus we arrive at a very important conclusion that, in the chiral limit, the mass-gap 
equation is the hound-state equation for the Goldstone boson at rest, responsible for the 
spontaneous breaking of chiral symmetry. Since the chiral pion is this Goldstone boson, 
then we are forced to see it already at the BCS level and it will remain being the massless 
Goldstone boson even when we go beyond BGS. Theproperties of the pion have been widely 
discussed in the literature — see, for example, for QCD in 3+1 and for QCD 

in 1+1. 

Beyond the chiral limit the solution to the bound-state equation takes an approxi- 
mate form: 



sm Lfp 



M^Ap 



Ml 



p^dp 
(27r)3 



Ap sin (y9p. 



(54) 



where all corrections of higher order in the pion mass are nedected and the function Ap 



obeys a reduced M7r-independent equation (see, for example, |^ or the papers [l^ where 
such an equation for Ap is discussed in two-dimensional QCD). 

Notice that, although the dressed quark dispersive law and T-amplitudes for the pion 
are real (see Appendix A for the details), the bound-state equation (jMj) admits a tachyonic 
solution with M'^^^^y < (see also the papers where a tachyonic solution was found for 
the bound-state equation in the trivial vacuum). For Mtach = "^l^tachl, the tachyon wave 
functions fail to be real and acquire imaginary parts, in accordance with Eq. It is easy 



to verify that they obey a simple transformation rule, (v^tach)* V^tach- 



3. The case of the p-meson 



Compared to the pion, the case of the p-meson involves an extra complication, a con- 
sequence of its richer spatial structure. The Hamiltonian ()25|1 does not commute with the 
operators of the total spin S and the angular momentum L, so it cannot be diagonalised in 
terms of the "^^^^Lj states. The appropriate basis is given by the set of physical observable 
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mesons J^^ . Then, although the lowest 0""*" state — the chiral pion — is represented by just 
only one term, a pure ^5*0 state, higher states are in general a mixture of angular momentum 
states, like, for instance, the case of 1 p-meson, which is a mixture of ^5*1 and ^Di terms. 
This problem does not appear in two-dimensional QCD and is a consequence of a richer 
structure of the spatial dimensions in the four- dimensional theory. Therefore, in order to 
formulate the eigenvalue problem for the p-meson, one should introduce a set of four wave 
functions, ^^tg^ = V^o ' V^sdi = ^^2 1' ^^"^ write the system of four coupled equation, 

q^dq 



[2E,-M,]ipt{v) 



{2.f 
q dq 

'W 

q dq 

w 

q dq 



+ ^02^ ^^^(9) + ^OO^V^O (?) + Tq2'V2{(1)] 

[Tti'^tiQ) + T^a^^o (?) + T^2'^2{q) + ^2V^o (?)] 

[T22^Vt{(l) + ^2"0^<^^(?) + ^2"2"^2 (?) + ^20"</'2 



(55) 



where the T-amplitudes can be built explicitly using the graphical rules established above 
and with the help of the spin-angular wave functions [K3s^{f))\siS2 = ['«o(^]siS2 
[k3£)^(^[s-^s2= [/t2(^]siS2- For example, similarly to Eq. we have. 



n2^{v.q) 



dVtpdVtqVip- q)Sp \ Ko{p)[u{p)Tu{q)]K2{q)[v{-q)Tv{-p)] 



(56) 



and so on. 

As an example, we give here the explicit form of the Bethe-Salpeter equation ()55p for 
the harmonic oscillator potential ^: 



1 





6 

+ - 


1 





'2 

+^ 





1 










2 



3-(Ti -2V2 
-2V2 3 + ai 



cos (pp 



2(1 -OTi) V2{l-a,) 
72(1 - ai) 1 - (Ti 



(T3 

-as 



2(1 — sin ipp) 

^2iP) 

.^oip) . 



-4 ^/2 
V2 -4 



(57) 



where ex's are the Pauli matrices and, similarly to the case of the pion, the new radial wave 
functions are defined as t'^2(p) — PVo2{p)- 

Naively, one may conclude, from the bound-state equation that the p-meson has to 
be described by a four-component wave function, rather than by just a two-component wave 



19 



function, as it was the case with the pion. This conclusion is erroneous, since the doubling 
of the number of equations in is a consequence of the usage of an inappropriate basis 
"^^^^Lj. Indeed, together with the light p-meson, we must have a second solution for the 
the bound-state equation needed to describe the heavier vectorial partner of the p. 
The easiest way to show this is to neglect off-diagonal terms in Eq. ()55|1 . thus splitting the 
system of four equation into two independent systems of two equations each, for yjg and iff, 
respectively. The eigenvalues found for these two independent eigenvalue problems give the 
masses of the pure ^Si and ^Di states, which should be mixed then as 

( M,s, - Mp AsD 

det 

y Ads Msd, - Mp 

with AsD and A^i^ being the contributions of the restored off-diagonal terms of the full 
Hamiltonian ()25|1 . The lighter solution represents the physical p-meson in the basis J^'^' . In 
other words, if an appropriate basis — which diagonalises the three-dimensional spatial part 
of the Hamiltonian — is chosen from the very beginning, each physical mesonic state is 
fully described by just a two-component wave function and the corresponding eigenvalue — 
the physical mass of the meson. Notice that the transformation ()58|1 can be viewed as yet 
another Bogoliubov-Valatin transformation parameterised by the S-wave-D-wave mixing 
angle. 

Now, with the p-meson included, the diagonalised Hamiltonian takes an obvious form, 

H = M^m|,m^ + MpUilmp + Mpiiv}p,mpi . . . , (59) 

where p' denotes the heavier partner of p which is also a solution of Eq. ()58|1 and the operators 
creating/annihilating the vector mesons are defined similarly to the pionic case, Eqs. ()27|) 
and ((221) 5 through the 1 wave functions '^^ pi- We are now ready to consider the general 
case. 




4. The general case 

Now we return to the general case of the Hamiltonian ()25j) and diagonalise it in terms 
of the mesonic compound operators. We work in the basis {n, J^*"}; n being the radial 
quantum number. We also assume that this basis includes all states with the given 
so that the set of such states is complete and has a one-to-one correspondence with the 
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complete set "^^^^Lj. For the sake of simplicity, we denote each mesonic state as {n, z/}, 
where the set of quantum numbers v is not only based on the J^*" classification scheme, 
but also identifies each meson in a given nJ^^ multiplet, as it was discussed in the case 
of p and p'. Then, for a given pair {n, i/}, representing a physical meson, we introduce a 
two-component radial wave function ^Pnuip) ^^"^ ^^e spin- angular wave function n^ip), as it 
was done in the pion case — see Eq. (j2Z|)- The operators and M of Eq.(j2nj, entering 
the Hamiltonian ()25|) . can be expanded using this basis (Eq. for the pion follows from 
this relation if we take only the lowest u = 0^^), 

Mss'{p,p) = 5^[/€.(^].s'M.(p). (60) 

V 

Now we use the completeness property of the set \^Ky{(p)] together with the fact that this 
set diagonalises the spin and the angular structure of the Hamiltonian ()25j) . to find: 

^ y d%\KXv)Xs^s-X>^v'(v)\s^s^= ^vv', (61) 



S\S2 



and 



^ / ci(]p(i(]g\/(p- g)[w++(p,g)],,,3,4,jK^(^]*2,jK^/(^],4,3 = -T++(p,g)5, 

S2S3S4 

/ dnpdn^V{p- q)[v~+{p,q)]s,s3SiS2[f^uii^]s^s2[fiu'CQ)]siS3 = -T~^{p,q)5. 

S2S3S4 

^ j d^pdngVip- g)[t;""(^,g)],,,3,4,jK^(^],,,2[K^.(^]*3,^ = -T--{p,q)5, 



(62) 



S1S2S3S4 



As a result, the Hamiltonian ()25|1 takes the form, in terms of the radial operators {Ml{p)} 
and {M^{p)} as: 

+Tt{P.<lWl{q)MUp) + T-+{p.q)MMM,{q) + r„-(p, ?)Mj(?)M„(p)} , (63) 

where all radial excitations are now disentangled from one another. 
Then the generalisation of the relations (fHUjl is trivial, 

I M^{p) = ^^[rfln^^nAp) + Hiu^nAv)] ^^^^ 
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f p'^dp ' 


j rhnu = J 


(27r)3 [ 
f p^dp ' 






(27r)3 [ 



Ml{p)^t{p) - MMVuAp) 

As in the pion case, we build the commutators of the operators rhnu and 

p^dp 



(65) 



p'^dp 



[ptu{p)vLu{p) - VnuiP)Vr, 



['^iuiP)'^nuip) - '^muip)'^LiP)] 



(66) 



(27r)3 

and require that they should obey the standard bosonic algebra, that is, [rhnu, ^Inu] = ^mn, 
and [rhnu, ^mu] = 0. Therefore, similarly to Eq. (jH^I) . for any given u, we arrive at the 
normalisation and orthogonality conditions for ip^s, ^ 

p'^dp 



(27r)3 
p^dp 



[ftuiPWtuiP) - fnuiP)^muiP)] = ^nm, 

(67) 

- '{'nui.P)'^'Lui.P)\ = 0- 



(27r)3 

The representation (jMj) . together with the normalisation and orthogonality conditions 
(|U7|). give the fully diagonalised Hamiltonian, 

n = Y^ Mnumi^mnu + O (-7^) , (68) 

n,u 

provided the mesonic wave functions obey the bound-state equation, 

[2E,-Mnu]vtu{p)= I ^[Tu^^ip,Q)^Liq) + Ttip,q)Vnuiq)] 

^^2! (69) 



q dq 
(2^ 



[2E, + Mnu]^nu{p) = I i^[T;^{p, q)y^L{q) + T;~{p, g)¥.;,(g)]. 



On comparing the mass-gap-like equation for the generalised Bogoliubov transformation 
with the Bethe-Salpeter equation pUj) . one can easily see that they coincide, the only 



"^Strictly speaking, the orthogonality condition for the mesonic wave functions should also include their 
spin-angular part. For example, 

^ j dnp[ti^Q)\l^,^[K^>$)\s^s^ j J^^['ptu{pW'Lu-iP) - 'Pnu{p)'P'^y'{p)\ =Snm.Suu', 

where, for non-coinciding i^'s, the l.h.s. vanishes due to the angular integration — see the orthogonality 
condition for k's, Eq. I|61fl . Then, for v — v' , one readily arrives at the first equation in ()67|l . 
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difference being that they are written in different representations: Eq. ()46p is written in the 
quark spins representation, whereas Eq. ()69|) is written in the J^^ representation. Therefore, 
it is a matter of rewriting the Salpeter wave function $^ in the proper representation to go 
from one equation to the other. Notice that the Lorentz nature of the confining interaction 
— the exphcit form of the matrices F in Eqs. and 1)411) — plays no role, defining only 
the explicit form the amplitudes v^^. Similarly, the Hamiltonian approach to the theory 
can be developed, following the same lines, for an arbitrary matrix F. 

In the leading order in A^^^; the Hamiltonian ()UH|) describes free stable mesons. The 
suppressed terms in flfjHj) must involve quark exchange and correspond to the parts of the 
Hamiltonian responsible for hadronic decays and scattering. For example, to consider a 
decay A B + C, one should restore the first sub-leading term in the Hamiltonian, 
-^=mA TTi^B with the coefficient giving the amplitude of the corresponding process 
12j |. Notice, however, that all three mesons — A, B, and C — cannot be at rest and one 



encounters a problem of boosting mesonic creation/annihilation operators, which is closely 
related to the general problem of Lorentz boosts in such potential models. This important 
issue lies beyond the scope of the present paper and deserves special treatment. Notwith- 
standing this general problem, we are still left with an important set of sub-leading physical 
processes which are amenable to exact treatment in the context of non-local NJL Hamil- 
tonians: the set of elastic hadron-hadron scattering at rest (scattering lengths) where we 
can define a common vacuum for all intervening hadrons [7j. We leave the issue of vacuum 
boosts to future publications. 

This completes the matching between the Hamiltonian and the Bethe-Salpeter ap- 
proaches to the non-local NJL models described in Eq.(^. We conclude that diagonalisation 
of the Hamiltonian of the theory in the mesonic sector in the leading order in Nc is equivalent 
to developing the Bethe-Salpeter approach to the theory in the ladder /rainbow approxima- 
tion. In the Hamiltonian approach, the mass-gap and the hound-state equations emerge to 
ensure the cancellation, in the Hamiltonian, of the anomalous Bogoliubov terms both in the 
quark and in the mesonic sectors of the theory. The normalisation condition for the mesonic 
wave functions ip^ with the minus sign between the positive- and the negative-energy com- 
ponents follows naturally in the Hamiltonian approach as the standard normalisation of the 
bosonic Bogoliubov amplitudes. 
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C. Discussion 

Let us make several concluding remarks concerning the diagonalisation procedure pre- 
sented in this section. We started from the Hamiltonian of a quark model with the in- 
stantaneous interaction parameterised by the quark kernel of an arbitrary, but necessarily 
confining, form. As an example for such like models we have, for instance, QCD in the 
truncated Coulomb gauge. After inclusion of self-interaction into the quark fields we obtain 
dressed quarks — the so-called BCS level. Chiral symmetry is spontaneously broken at this 
stage, quarks acquire effective mass and the chiral condensate appears in the vacuum. The 
next step was to perform yet another Bogoliubov-type transformation, so as to build oper- 
ators which create/annihilate nJ^^ mesonic states in the BCS vacuum. Finally, the rules 
relating the second Bogoliubov transformation (which diagonalised, up to A'^p-suppressed 
terms, the quartic part of the Hamiltonian) and the Bethe-Salpeter equation for bound 
quark- ant iquark states were deduced. Crucial points of the approach necessary to carry out 
this program were: 

• instantaneous inter-quark interaction which allows one to avoid the problem of the 
relative inter-quark time and to formulate a self-consistent Hamiltonian approach to 
the theory; 

• limit of a large number of colours with the inherent suppression of all non-planar 
diagrams (quark exchange); 

• confinement, which allows one to reformulate the theory entirely in terms of colourless 
bound states — mesons (the crucial conjecture (j23|l will fail for nonconfining interac- 
tions); 

• a nontrivial solution to the mass-gap equation defining the broken phase of the the- 
ory. This phase is characterised by two concurrencial processes, one defined by the 
amplitude of creation of a quark-antiquark pair by means of the 1)^3) operator applied 
to the vacuum (M]^,(p, jf)), and a second process realised through the "borrowing" of 
qq pairs from the chiral condensate via the operator dh (Mss'iPyff))- Both processes 
are fundamentally important for the chiral pion to become massless due to strong can- 
cellations between the positive- and negative-energy components of the pionic wave 
function. 
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The discussed approach is completely insensitive to: 

• the number of spatial dimensions, provided a suitable basis is built which diagonalises 
the spatial part of the Hamiltonian. This is due to the fact that both the mechanism 
for the separation of positive from negative-energy components of the mesonic wave 
functions and the mechanism for the disentanglement of radial excitations are universal 
mechanisms for instantaneous interactions; 

• the Lorentz nature of confinement, which we considered to take the simplest form — 
7o X 7o, but which can be easily generalised to be 'jfj, x 'ju', 

• the form of the potential, which is required only to be confining and to lead to a 
finite mass-gap equation. All power-like potentials with ^ a ^ 2 meet these 
conditions. Among these we have the linearly rising potential, a = 1, — as the most 
natural candidate for confinement in QCD, as well as the harmonic oscillator potential, 
a = 2, which is the easiest example for analytical and numerical studies. 

Finally, we arrive at the following important conclusions: 

1. the Bethe-Salpeter equation for bound states of quarks and antiquarks is equivalent 
to diagonalisation of the quartic term in the quark Hamiltonian, expressed in terms 
of the dressed quark fields. Notice that no new information, at least in the leading 
order in Nc, is required/appears beyond BCS level, and then we naturally arrive at 
the second conclusion that 

2. the entire problem is completely defined as soon as the mass-gap equation is formulated 
and solved, the chiral angle being the only entity which one needs to know in order to 
solve the entire theory. In particular, 

3. the second, bosonic, Bogoliubov transformation is defined via the Salpeter amplitudes 
ip~^ and ip~ . In the literature of boson condensation they are usually denoted by u 
and V amplitudes respectively. Therefore it is not hard to understand that ip^ play 
simultaneously the role of the mesonic wave functions and satisfy the normalisation 
condition which is proper of the usual Bogoliubov self-consistency condition u^—v'^ = 1. 
Thence these two amplitudes can be parameterised as m = coshx and v = sinhx, 
with X being the mesonic Bogoliubov angle. Therefore, for the generalised mesonic 
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transformation, cosh x and sinh x are given by the sets of the positive- and negative- 
energy components {'■Pny} and {'{>nu\ which are completely defined by the chiral angle 
(for example, for the chiral pion we have = —ip:^ = J^f^sinipp). 

In short, solving the mesonic Bethe-Salpeter equation is tantamount to finding the bosonic 
Bogoliubov transformation to fully diagonalise — in the leading order in — any non-local 
NJL Hamiltonian. 

This ends the discussion of the diagonalisation of non-local NJL models. In the next 
section, we generahse the results of this section to the case of multiple vacuum states — 
replicas. 

III. VACUUM REPLICAS IN POTENTIAL QUARK MODELS 

A. Multiple solutions to the mass-gap equation 

As noticed above, the mass-gap equation is equivalent in momentum space to a 
Schrodinger-like equation, with the dressed quark dispersive law playing the role of the 
effective potential (see the example of Eq. (fTH|) ). This suggests that the solution to the 
mass-gap equation might not be unique. Such multiple solutions were discovered for the 
harmonic oscillator potential in . Recently a detailed analysis was performed for the 
general form of the power-like confining potential, Eq. Q, and the existence of an infinite 



tower of solutions to the corresponding mass-gap equation was proved Let us enumerate 



the main statements one can make concerning the mass-gap equation: 

• in order to provide a nontrivial solution to the mass-gap equation the dressed quark 
dispersive law must become negative in the infrared region p 0. Since the chiral 
symmetry is broken only for such nontrivial chiral angles, then this property of Ep is 
absolutely necessary for chiral symmetry breaking; 

• once a nontrivial solution ipp to the mass-gap equation is found, it defines the wave 
function of the pion — the Goldstone boson — as ip = M-^smipp] 

• the slope of the chiral angle at the origin defines the scale of the chiral symmetry 
breaking for the given solution; steeper solutions correspond to less broken symmetry 
and to a sharper behaviour of Ep at the origin; 




26 



• the mass-gap equation for any single-parameter confining potential leads to an infinite 
number of solutions for the chiral angle. 

To exemplify the appearance of the replicas, let us draw the following qualitative picture. 
For the trivial chiral angle ipp = there is no dressing of quarks, the quark dispersive law 
being just Ep = p for all momenta. With such Ep, the effective potential in the Schrodinger 
equation (and similarly for other forms of the confining potential), Vesip) = '2Ep{p), 
does not possess a single bound state with a zero eigenvalue — therefore, no nontrivial 
solutions to the mass-gap equation do exist. We encounter here the well-known problem 
of the constituent quark models. Indeed, as it was discussed above, the mass-gap equation 
plays the role of the bound-state equation for the pion at rest. Upon substituting the free 
quark dispersive law, Ep = p, into the mass-gap equation — see for example, Eq. (fT5|) — 
we are led to the usual Schrodinger equation for the pion, [2p + V{r)]ilj.„ = M-^ip.,,, which 
is characteristic of potential constituent quark models. The pion mass coming out of this 
equation appears to be of order of the confining interaction scale, ~ 400 -^ 500Me^, 
which is several times larger than the experimental value of lAOMeV and such M^. does not 
vanish in the chiral limit. In the bound-state equation, this problem is known to be solved 
by the presence of the negative-energy component ip~ of the pionic wave function which 
happens to be of the same order of magnitude as the positive-energy component Strong 
cancellations between the two components of the pion wave function bring the pion mass to 
zero. At BCS level, the solution of this problem comes from the actual form of the dressed 
quarks dispersive law Ep which becomes negative at small momenta, and cancels the positive 
contribution of the confining potential allowing for the pion mass to vanish. In other words: 
the peculiar behaviour of the dressed quark dispersive law is both a necessity and a direct 
consequence of the chiral symmetry breaking. Suppose now that we can parameterise this 
small-p negative contribution to Ep by a mass scale /i. For the power-like potentials (jH)), it is 
clear that such a contribution has to be proportional to Kq'^°' (see the definition of Ep through 
the auxiliary functions Ap and Bp, Eq. ()12j). as well as Eqs. and (fTUI) ). For dimensional 
reasons, Ep{p = 0) = —const °^ (for example, for a = 2, from Eq. (fT^ with m = and 
(Pp ~^ f — ^ one readily finds that Ep{p = 0) = — f-^)- Thus, for a sufficiently small /i = /io 
the effective potential V^ff(p|/i) becomes binding enough to produce an eigenstate with a zero 
eigenvalue. Then /xq defines the scale of the chiral symmetry breaking in the new vacuum 



27 



and the corresponding chiral angle behaves as !fo{p) ~ ^ — — + . . .. This is the ground-state 
solution which corresponds to the BCS vacuum of the theory. Decreasing the scale fi, one 
can reach a situation where a second zero eigenvalue bound state appears for the potential 
Ves{p\fJ'). According to general quantum mechanical theorems, this solution must contain 
one knot. This is the first replica vacuum. Continuing this procedure, we can build the 
whole infinite tower of replicas for a given potential. In other words, for any given value of 
the scale /i one has a set of orthogonal eigenstates (for the potential Kff(p|/u)) with positive 
and negative eigenvalues. To search for the nth replica, it is sufficient to vary fi, thereby 
shifting the whole tower of eigenstates up or down until the appearance of the sought nth 
zero eigenvalue. Therefore, different replicas are eigenstates in different potentials. 

The easiest way to prove such a picture is to evaluate the quasi-classical Bohr-Sommerfeld 
integral for such an eigenvalue problem. Since we consider a one-scale confining interaction, 
then, from dimensional reasons, it is clear that the kinetic and the potential energies, in mo- 
mentum space, behave as J^i+"r° and Kq'^" /p", respectively. Consequently the correspond- 
ing WKB inte gral depends logarithmically on the scale n, which plays the role of the cut-off, 
-^WKB oc In Therefore, the quasi-classical quantisation condition, Iwkb = 27rn, can 
be fulfilled for any n, provided the corresponding scale /i„ is small enough. For example, for 
the harmonic oscillator potential in D dimensions, the approximate dependence of the scale 
fi on the index of the replica n can be easily found using the aforementioned WKB method 
to be. 



so that, in any dimension D > 2, the number of replicas is, indeed, infinite. The boundary 
case of D = 2 contains a singularity, for n 7^ 0, so there are no replica solutions in the 
harmonic oscillator potential in two dimensions, as was found in and discussed in detail 



Of course, the procedure of building replicas drawn above should be understood as a 
simplified qualitative method, since the scale fi is not a free parameter but it has to be 
defined self-consistently with the form of the chiral angle. Notice also that, for confining 
potentials defined through more than one parameter, the proof given above does not hold, 
since the WKB integral may be regularised by the second scale, instead of the parameter 
fi. Then the number of replicas becomes finite or they do not exist at all. Notice, however. 




(70) 



in ll7|. 
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that, in the general case, it is much harder not to have rephcas than to have them, and the 
expectation of existence of more than one solution for nonlinear mass-gap equations is not 
only physically attractive but also quite natural from the mathematical point of view. 

B. The problem of the tachyon 



In this subsection, we address the problem already touched upon in the paper 




namely, the problem of the tachyon which appears in the excited vacua. Indeed, it was 
numerically found that, in the chiral limit, the chiral condensate changes the sign from replica 
to replica remaining negative for even states (the ground-state BCS vacuum being the lowest 
representative of the even states) and becomes positive for all odd states, starting with the 
first replica 0, The universal status of this rule becomes clear from the asymptotic 
behaviour of the chiral angle (fTTj) and the fact that for each next replica the chiral angle 
possesses an extra knot and, therefore, it approaches, when p — oo, the (fp = asymptote 
from the opposite half-plane, as compared to the previous replica. Then, switching on a 
small quark mass, we arrive at the Gell-Mann-Oakes-Renner relation ()18|). in which the 
chiral condensate on the r.h.s. can be calculated in the limit m = ^, and, for positive 
values of (qq), the pion becomes a tachyon, < 0. As always, the presence of a tachyon 
means that the lowest state — the vacuum — is chosen improperly, and there should be a 
more preferable vacuum. It is easy to see that, for odd replicas, this is indeed the case. 

To demonstrate this, let us consider the exact chiral limit and the most general Valatin- 
Bogoliubov transformation from the trivial vacuum |0)o to the generalised vacuum \9) 
which can be written as, 

\e) = So\0)o, Se=exp[Ql-Qe], qI = ]^Y.Y. I 4^^^^"^^^'^^^^^^^ (^1) 

a ss' ^ ' 

where (/jp is the chiral angle and 

m,e\ = m^p^ cos 9 + miso sin 9. (72) 



■^Beyond the chiral hmit the quark condensate acquires an infinite contribution coming from the trace 
of the free-particle Green's function which should be subtracted. The term of the zeroth order in m in the 
properly regularised condensate obviously coincides with the value calculated in the chiral limit. 
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The ^Pq matrix Tlap^ = {af))ia2 is studied in detail in ^ and the ^5*0 matrix 97li5g = 
DJItt = —io'2- For = Eq. (|7T|l reproduces the standard definition of the BCS vacuum |0) 



. When 6* = I we obtain: 



2 / = 5*^/210)0 = exp 



|0)o = exp[zg5]|0)o, (73) 



with the operator being responsible for translations along the 6 direction. The angles 
9 and ipp are independent quantities. For a given solution of the mass-gap equation, ipp, 
when = TT we get S',r[v^p] = So[—ipp], which, as discussed above, is the usual pseudo-unitary 
^Pq operator for the chiral condensation. Thus, translation from ^ = to ^ = vr along the 
6 direction is tantamount to a chiral angle transformation of (pp —fp- We could have 
started with —ipp — also a solution to the mass-gap solution for rrig = — to arrive, in the 
end of the ^-journey, at (fp. Notice that, in general, the state |^) does not have any definite 
parity — we have to consider mixtures of different ^-vacua in order to build a state with a 
given definite parity. For instance, \G) + \ — d) has the parity plus. In the special case of 
6 = 71, = S-TT and the state = tt), as well as = 0), has positive parity. 

Let us start with a given ipp and put, for a moment, 6 = 0. If we now plot the vacuum 
energy as a function of the chiral condensate S = (qq), then S = corresponds to the local 
maximum, whereas the stable minimum is provided by a nonzero value of S, easily related 
to the corresponding chiral angle, 

Nc 



7r2 



dp p sin ipp. (74) 



Next, we can move along the 0-valley. Since, in the chiral limit, the charge commutes 
with the Hamiltonian of the theory, [Q5H] = 0, then the vacuum energy is degenerate for 
all ^'s, so that we are free to choose any value for 6. Such a form of the vacuum energy as a 
function of two variables, S and 6, is known as the Mexican hat. Incidentally, we can find, 
starting from expression (fTHj) . the chiral pion Salpeter amplitude (p^{p) = —(p^{p), which is 
just the c-number multiplying the anomalous part of ,18]- It is then clear that the pion 
Salpeter amplitude is given by simpp, independently of any consideration for a given quark 
kernel, provided it supports the mechanism of spontaneous chiral symmetry breaking. 

If now a small quark mass is introduced, then the vacuum energy acquires a contribution 
of the chirally non-invariant term proportional to the quark mass m and, as a result, the 
whole picture gets tilted — the state with the negative sign of the chiral condensate being 
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FIG. 2: The ground-state (the left plot) and the first replica (the right plot) solutions of the mass- 
gap equation after the redefinition of the chiral angle sign for odd replicas. The current quark 
mass is m = O.Ol-f^o. The momentum p is given in the units of Kq. 

energetically preferable as compared to the state with the positive condensate. These two 
states differ by their 6 coordinate — one of them corresponds to 6 = 0, whereas the other 
one has 6 = n. The choice of the chiral angle adopted before, with ipp{0) = ^, selects the 
vacuum 6' = as the true vacuum of the theory |0). This is indeed the case for the ground 
BCS vacuum, as well as for all even replicas, with S < 0. The pion, which is not tachyon in 
these vacua, is responsible for the vacuum energy increase during the rotation in the angle 
6 from |0) to |7r). For odd replicas, the chiral condensate changes the sign, so that the state 
with 6 = IT must acquire a lower energy than the one with ^ = 0. As a result, in odd replicas, 
the pion becomes the tachyon responsible for the energy decrease during the same rotation 
form |0) to \tt), the latter representing now the true vacuum. Therefore, for odd replicas, 
we set V9p(0) = — | instead of |. This choice will ensure that the chiral angle, when the 
momentum p goes to infinity, will approach zero from above. Then the chiral condensate 
flTIj) changes the sign, and so does the pion mass squared — the pion mass becomes real. As 
discussed before, the definition of the chiral angle admits such a change and, as a result, we 
arrive at two different classes of solutions to the mass-gap equation: even solutions, which 
start from | at p = and approach the free limit at large p's from above, and odd solutions 
which also approach their large-p asymptote from above but, at p = 0, start from — ^. This 
solves the problem of the tachyon for odd replicas, making the latter normal vacuum states 
with the possibility of building the spectrum of hadrons above them. 

As an example, we give in Fig. 2, the profiles of the ground-state and the first replica 
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FIG. 3: Solutions of the bound-state problem for the pion, Eq. (|?T|) . in the BCS vacuum (the left 
plot) and in the first replica with the reversed sign of the chiral angle (the right plot). The overall 
normalisation factor is omitted, the largest component being fixed equal to unity at p = 0. The 
current quark mass is m = O.OIKq. The momentum p is given in the units of Kq. 

solutions to the mass-gap equation (fT^ for the oscillator-type potential (the case of the 
generalised power-like confining potential ^ is studied in detail in [l^)- discussed 
before, the sign of the chiral angle for the first replica vacuum is reversed. 

In Fig. 3, we represent the solutions of the bound-state problem for the pion with the 
harmonic oscillator confining potential, Eq. (j51|) . for the two lowest solutions to the mass- 
gap equation (fT^ depicted in Fig. 2 — for the ground BCS vacuum (the left plot in Fig. 3) 
and for the first replica with the reversed sign of the chiral angle (the right plot in Fig. 3). 
It is clearly seen from Fig. 3 that the pionic wave functions are indeed very close to ± sin ipp 
(see Fig. 2) with the corrections given by the formula (jSlj) . Notice that, after the change 
(fp —(fp performed for odd replicas, the "+" and the "— " components of the mesonic 
wave function substitute each other, so that one needs either to redefine the whole set of the 
operator transformations used to diagonalise the Hamiltonian of the theory in the mesonic 
sector or, which is more economical, simply to rename ip^ to (p~ and vice versa. 



C. Hadronic spectrum in the replica vacua 

It is concluded, in the second section, that any NJL-type Hamiltonian model with an 
arbitrary confining quark kernel is completely defined at BCS level, as soon as the chiral 
angle is built. If the large- A^^c limit is assumed, then all approximations become controllable 
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and the Hamiltonian of the theory admits complete diagonahsation in terms of the compound 
mesonic operators. We found that the appropriate basis for such a diagonahsation was 
provided by the set {n, J^^} and each mesonic state was described by a pair of wave 
functions y?^. In this section, following the papers P, l^ll^, we gave evidence of existence 
of excited solutions to the mass-gap equation — the vacuum replicas. Now we address the 
question concerning the spectrum of hadrons in the replica vacua. 

As argued in the previous subsection, the pion is massless in the chiral limit for any 
replica vacuum and, with the properly defined chiral angle, it acquires a small real mass 
beyond the chiral limit. As far as highly excited mesonic states are concerned, the chiral 
symmetry is restored in this part of the spectrum and the hadrons' properties become 
insensitive to the details of the vacuum in which they are created since, for high excitations, 
the negative-energy components ^Pnuip) wave functions are negligible and the Bethe- 

Salpeter bound-state equation (jU^ amounts to a quantum-mechanical Schrodinger equation. 



Consequently, no cancellations between the positive- and negative-energy components take 
place anymore, as it happens to the chiral pion and leads to the small pion mass of lAOMeV, 
as opposed the value of about AOO-rbOOMeV which would follow from Eq. (|75|). As a result, 
in the vacua with "less" broken chiral symmetry (higher replicas), the masses of mesons are 
pushed up, with the exception of the chiral pion, whose mass, in the chiral limit, is kept 
equal to zero by the requirement of the chiral symmetry. Beyond the chiral limit, the pion 
mass also increases with the index of the replica (see Fig. 3 and compare the pion mass of 
0.57Kq, in the BCS vacuum, with the value of 0.988/^0; for the first replica). 

Therefore, we conclude that each mesonic state can be characterised by an extra "quan- 
tum number" — the index of the vacuum replica in which it was created. Then the full 
diagonalised Hamiltonian of the theory takes the form: 



Af being the replica label, which is the generalisation of Eq. ()68|) for the multi-vacuum case. 





(75) 




(76) 
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IV. CONCLUSIONS 

In this paper, we have shown how the generahsed NJL-type Hamiltonian model with a 
nonlocal confining quark kernel can be fully diagonalised in the sector of physically observable 
hadronic states. We extend this result to include the multiple vacuum states recently found 
to exist in models of such a type and which are very likely to exist in real QCD. The main 
result of the paper is the fully diagonalised Hamiltonian ()76|) . Similarly to the case of a path 
integral with a multi- minimum function in the exponent, the Hamiltonian (fTBj) sums the 
contributions of all minima, with the corresponding weight — the masses of mesons in the 
given vacuum — being minimal for the lowest vacuum state and increasing for higher vacua, 
thus suppressing the contribution of the corresponding replicas to the full partition function 
of the theory. We study the form of the chiral angle profiles defining the excited vacua 
and find that for an even (odd) replica to decay to a lower even (odd) one it is sufficient 
to change the chiral angle in the finite region in momentum p, since both asymptotes for 
both replicas of the same parity coincide: ^Pp{p = 0) = | (^f); Vp{P ^ = 0- the 
contrary, for the decay of an odd replica to an even one, or vice versa, the chiral angle has to 
be changed dramatically in the infrared region, since the low-momentum behaviour of even 
and odd replicas is quite different. 

In the paper [9|] a quantum field theory perturbative method was developed. It used 
the replica-quark-antiquark vertex parameterised by a small difference between the chiral 
angles of the ground-state vacuum and the rephca to evaluate the quark self-energy shift 
due to the replica presence. This method can be easily extended to the case of an infinite 
number of replicas, but notice that the sets of even and odd states have to be considered 
separately. As discussed above, only transformations of the chiral angle between states with 
the same parity are small and can be considered perturbatively. We conclude, therefore, 
that two sets of perturbative approaches to replicas should be built — for even and for 
odd replicas, independently. On the contrary, a transition between states with opposite 
parities, for example, the decay of the first replica to the ground BCS vacuum, involves 
the global transformation of the chiral angle generated by the pseudoscalar pionic operator 

and results in a burst of pions with a huge energy release. The same transition in the 
opposite direction — excitation of the replica — requires an external global source which is 
not inherent to the model. Building of such a source is an important problem in the theory 
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of replicas and it will be the subject of future publications. 

Acknowledgment s 

One of the authors, A. Nefediev, is grateful to P. Bicudo and Yu. S. Kalashnikova for many 
fruitful discussions and would like to thank the staff of the Centro de Fisica das Interacgoes 
Fundamentals (CFIF-IST) for cordial hospitality during his stay in Lisbon, where this work 
was originated, and to acknowledge the financial support of the grant NS-1774.2003.2, as well 
as of the Federal Programme of the Russian Ministry of Industry, Science and Technology 
No 40.052.1.1.1112. 



[1] Y. Nambu, G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 

[2] A. Amer, A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. Raynal, Phys. Rev. Lett. 50, 87 (1983); 

A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. Raynal, Phys. Lett. 134B, 249 (1984); Phys. 

Rev. D 29, 1233 (1984). 
[3] A. Le Yaouanc, L. Oliver, S. Ono, O. Pene and J.-C. Raynal, Phys. Rev. D 31, 137 (1985). 
[4] P. Bicudo and J. E. Ribeiro, Phys. Rev. D 42, 1611 (1990); ibid., 1625 (1990); ibid., 1635 

(1990); P. Bicudo, Phys. Rev. Lett. 72, 1600 (1994); P. Bicudo, Phys. Rev. C 60, 035209 

(1999) . 

[5] N. H. Christ and T. D. Lcc, Phys. Rev. D 22, 939 (1980); see also A. P. Szczepaniak, E. S. 
Swanson, Phys. Rev. D 65, 025012 (2002) and references therein. 

[6] S. L. Adler and A. C. Davis, Nucl. Phys. 244B, 469 (1984); Y. L. Kalinovsky, L. Kaschluhn, 
and V. N. Pervushin, Phys. Lett. 231B, 288 (1989); P. Bicudo, J. E. Ribeiro, and J. Rodrigues, 
Phys. Rev. C 52, 2144 (1995); R. Horvat, D. Kekez, D. Palle, and D. Klabucar, Z. Phys. C 68, 
303 (1995); Yu. A. Simonov, Yad. Fiz. 60, 2252 (1997) [Phys. Atom. Nucl. 60, 2069 (1997)]; 
N. Brambilla and A. Vairo, Phys. Lett. 407B, 167 (1997); Yu. A. Simonov and J. A. Tjon, 
Phys. Rev. D 62, 014501 (2000); P. Bicudo, N. Brambilla, E. Ribeiro, and A. Vairo, Phys. 
Lett. 442B, 349 (1998); F. J. Llanes-Estrada and S. R. Cotanch, Phys. Rev. Lett. 84, 1102 

(2000) . 



35 



[7] P. Bicudo, S. Cotanch, F. Llanes-Estrada, P. Maris, E. Ribeiro, and A. Szczepaniak, Phys. 
Rev. D 65, 076008 (2002). 

[8] P. J. A. Bicudo, A. V. Nefediev, and J. E. F. T. Ribeiro, Phys. Rev. D 65, 085026 (2002). 

[9] A. V. Nefediev and J. E. F. T. Ribeiro, Phys. Rev. D 67, 034028 (2003). 
[10] P. J. A. Bicudo and A. V. Nefediev, Phys. Rev. D 68, 065021 (2003). 
[11] A. A. Osipov and B. Hiller, Phys. Lett. 539B, 76 (2002). 

[12] Yu. S. Kalashnikova and A. V. Nefediev, Yad. Fiz. 62, 359 (1999) [Phys. Atom. Nucl. 62, 323 
(1999)]; Usp. Fiz. Nauk 172, 378 (2002) [Phys. Usp. 45, 347 (2002)]; Yu. S. Kalashnikova, A. 
V. Nefediev, A. V. Volodin, Yad. Fiz. 63, 1710 (2000) [Phys. Atom. Nucl. 63, 1623 (2000)]. 

[13] G. 't Hooft, Nucl. Phys. B75 (1974) 461. 

[14] K. Kikkawa, Ann. Phys. 66, 3633 (1981); A. Nakamura and K. Odaka, Phys. Lett. 105B, 
392 (1981); Nucl. Phys. 202B, 457 (1982); S. G. Rajeev, Int. Journ. Mod. Phys. A9, 5583 
(1994);A. Dhar, C. Mandal, and S. R. Wadia, Phys. Lett. 329B, 15 (1994); A. Dhar et.al. Int. 
Journ. Mod. Phys. AlO, 15 (1995); M. Cavicchi, Int. Journ. Mod. Phys. AlO, 167 (1995); K. 
Itakura, Phys. Rev. D 54, 2853 (1996). 

[15] Yu. A. Simonov, Yad. Fiz. 60, 2252 (1997) [Phys. At.Nucl. 60, 2069 (1997)]; N. Brambilla 
and A. Vairo, Phys. Lett. 407B, 167 (1997); Yu. S. Kalashnikova and A. V. Nefediev, Phys. 
Lett. 414B, 149 (1997). 

[16] M. Gell-Mann, R. J. Oakes, and B. Renner, Phys. Rev. 175, 2195 (1968). 

[17] P. J. A. Bicudo and A. V. Nefediev, Phys. Lett. 573B, 131 (2003). 

[18] P. J. A. Bicudo, J. R. Rodrigues and J. E. F. T. Ribeiro, Phys. Rev C 52, 2144 (1995). 

[19] I. Bars, M. B. Green, Phys. Rev. D 17, 537 (1978). 



APPENDIX A: THE PION BETHE-SALPETER EQUATION 

In this appendix we derive the Bethe-Salpeter equation for the pion for the generic form 
of the potential and for the Lorentz structure of the confinement being 70 x 70. We generalise 
the method suggested in the paper ^ for two-dimensional QCD. We start from Eq. (j4ip 
for the mesonic Salpeter amplitude and define a matrix mesonic wave function, 

^{p; M^) = j ^S{p,po + Mj2)x{p- M^)S{p,po - Mj2). (Al) 
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We also present the Dirac projectors ()37p in the form: 



(A2) 



and introduce a modified wave function ^{p;M^) = T^^ {p; M^)T^ . Equat ion for this new 
matrix wave function following from Eq. ()41|) . with F = 70, reads: 



^{p;M^) 



V{p- q) 



^ TlT,^{q;M^)T,Tl ^ ^ ^ T^T,^ jq; M^)T,T^ 



2Ep - 



2Ep + 



It is clear that a solution of Eq. ()A3|1 has the form, 



(A3) 



(A4) 



and, due to the obvious orthogonahty property of the projectors P±, P+P- = P-P+ = 0, only 
matrices anti-commuting with the matrix 70 contribute to A and B. The set of such matrices 
is {75,7075,7,707} which can be reduced even more, up to {75,7}, since the matrix 70 can 
be always absorbed into projectors P±. For the case of the chiral pion only 75 contributes 
and one has: 

= Ib^tip), = -l5Vn{p), (A5) 

where the signs and the coefficients are chosen such that to comply with the definitions ()29|l 
and (jini)- It is an easy task now to extract the amplitudes (see Eq. ((211)) from Eq. ()A3|) 
using the explicit form of '^ni'P', M^,) and the operator Tp. They read: 



dn,V{p~ q) 
dn,V{p- q) 



COS' 



sm 



COS COS 



(Pp-iPq 1 - 



+ 



COS (Pp COS (Pq 



(A6) 

For the harmonic oscillator potential, and with the amplitudes ()A6|) . Eq. (|34j) reproduces 
Eq. (jni)- 

In the chiral limit, f^{p) = —ip~{p) = v^7r(p), so that the bound-state equation 
reduces to a single equation. 



2Epip^{p) 



q^dq 
(2^ 



[T++(p,g)-T+-(p,g)]^.(g) 



d?q 
(2^ 



V{p-q)^M. (A7) 
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or, in the coordinate space, one arrives at the Schrodinger-hke equation. 



[2Ep + V{r)]^^ = 0. 



(A8) 



It is also instructive to derive the bound-state equation for the pionic matrix wave function 
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\l/(p;M7r). It follows directly from the representation ()A3|) after a simple algebra (see 
for the detailed discussion of the matrix wave function formalism in two-dimensional QCD), 



M,^(p; M,) = [{dp) + 7om]^(p; M,) + M^)[{ap) - 70m] 



+ 



(A9) 



-A-(g-)^(p;M0A+(-g) + A-(p)*(g1M^)A+(-p)}. 

The explicit form of M^r) through the components ip^ and the chiral angle can be 
found easily from Eqs. ()A4|1 . ()A5j) . 



(AlO) 



where = i(<^+ - ), = i(<^+ + ). 

Let us multiply the matrix bound-state equation ()A9|) by 7075, integrate its both parts 
over d^p, and, finally, take the trace. The resulting equation reads: 



d^p 



sin = 2m 



d^p 



(All) 



It is easy to recognise the Gell-Mann-Oakes-Renner relation ()18|) in Eq. ()A11|) . Indeed, 
using the explicit form of the pionic wave function beyond the chiral limit, Eq. ()54|) . one 
can see that = (AC/ V M^) sin ipp and = (ACv^^^^) Ap which, after substitution to 
Eq. (jAlljl . give the sought relation. 



Ml I dpp'^ApSiii.Lfp =-2m{qq), 

Jo J 

where the definition of the chiral condensate (j74p was used. 



(A12) 



